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Abstract 

We observe that a recently proposed supersymmetric model with flavor symme- 
try admits a new CP violating ground state. A new sum rule for the quark mixing 
parameters emerges, which is found to be consistent with data. Simple extensions of the 
model to the neutrino sector suggest an inverted hierarchical mass spectrum with nearly 
maximal CP violation (I^mnsI — '^l'^)- Besides reducing the number of parameters in 
the fermion sector, these models also provide solutions to the SUSY flavor problem and 
the SUSY CP problem. We construct a renormalizable scalar potential that leads to 
the spontaneous breaking of CP symmetry and the family symmetry. 
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1 Introduction 



Non-Abelian discrete symmetries have found applications in explaining aspects of the fla- 
vor question not addressed by the standard model (SM) of particle physics. Restrictions 
imposed by such symmetries can lead to predictions for the Cabibbo-Kobayashi-Maskawa 
(CKM) quark mixing angles in terms of the quark mass ratios [Ij. Such symmetries have 
been employed successfully to generate a geometric structure for the leptonic mixing angles, 
independent of the lepton mass ratios [2]. 

The non-Abelian discrete symmetric structure G appears to bode well within the su- 
persymmetric (SUSY) standard model, since the same symmetry can also provide a natural 
solution to the excessive flavor change that occurs with generic soft SUSY breaking terms [3]. 
The three families of quarks and leptons will transform as doublets or triplets of the group 
G, which would result in the degeneracy of their masses. Degenerate squarks (and sleptons) 
would alleviate the SUSY flavor violation problem. If the non-Abelian symmetry is con- 
tinuous [1] and gauged [3], there is flavor violation arising from the D-terms of the flavor 
group [H]. Continuous global symmetries are susceptible to explicit violation from quantum 
gravity. Non-Abelian discrete symmetries which have a gauge origin are free from these 
problems and deserve special considerations [7]. 

In Ref. [H] we presented a SUSY model based on the non-Abelian discrete group Qq - the 
binary dihedral group of order 12. This group has two inequivalent doublet representations, 
one real doublet, and one pseudo-real doublet, which can be handy for model building. In the 
flavor sector this symmetry results in one prediction for for a combination of the CKM mixing 
parameters, which was shown to be consistent with data. The quark and lepton superfields are 
assigned to doublets and singlets of Qq, with the singlets identified as belonging to the third 
family. In the symmetric limit the squarks of the first two families would be degenerate, 
which is sufficient to solve the SUSY flavor problem. Furthermore, by assuming that CP 
violation has a spontaneous origin, this model also solves the SUSY CP problem. Excessive 
CP violating processes arising from the SUSY breaking sector are absent, since the parameters 
are all real. Yet the model admits CP violation in the quark mixing matrix. 

One major purpose of the present paper is to show that the Qe model studied in Ref. [8] 
admits a new minimum which violates CP, but leaves a new interchange symmetry in tact. 
By virtue of this interchange symmetry, we derive a new sum rule among the quark mixing 
parameters and CP violation, which is found to be consistent with observations. Such an 
interchange symmetry was present in Ref. |8] as well, but the new one presented here is 
different, although it arises from the same Higgs potential. We extend this symmetry to the 
lepton sector and obtain interesting correlations between the neutrino oscillation parameters. 
We also compare the predictions of the new minimum with those of the old, and in the process 
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update our old predictions. We use the most recent values of light quark masses where the 
errors have decreased significantly as a result of improved lattice calculations. We compare 
the model predictions to the best fit values obtained in the SM as well as by including certain 
new physics contributions in Bd^s — Bd,s mixings as obtained by the CKFfitter group [9]. 
These new contributions are motivated by certain discrepancies obtained in the SM CKM 
fits - such as the differences of order 20% in the CP violation parameter rj obtained from 
fits to ex and B — > J /'^Ks decay. Small new physics contributions naturally arise in our Qq 
based model. For example, there are contributions to meson-antimeson mixing via SUSY box 
diagrams, which may be important for the Bd^s meson system since the third family squark 
is not degenerate with the first two family squarks. 

We also present a complete Higgs potential that leads to the spontaneous breaking of CP 
symmetry and the Qq fiavor symmetry without leading to pseudo-Nambu-Goldstone bosons. 
In addition to the Qg, a fiavor universal symmetry is introduced. Owing to this Z4, even 
after spontaneous symmetry breaking, an unbroken interchange symmetry survives in the 
Higgs potential, for which there are two possible choices, denoted as P/,//. These symmetries, 
along with Qe, reduce significantly the number of parameters in the fermion mass matrices. 
This reduction of parameters leads to a sum rule involving quark masses and mixings [8]. 
Moreover, CP violation has a spontaneous origin, which is perhaps more satisfying than the 
usual assumption of explicit CP violation. Nevertheless, the dominant source of CP violation 
in the quark sector is the Kobayashi-Maskawa mechanism. The sum rule involving quark 
masses and mixings that has been derived relies on the spontaneous violation of CP. With 
this, the problem of excessive CP violation that generically exists in the soft SUSY breaking 
sector can be solved in a rather simple way. Various phenomenological aspects of this model 
in minimum Vi have been studied in Ref. [Sl fTUlfTT] . 

The plan of the paper is as follows. In Sec. 2 we present the supersymmetric Qq model. 
In Sec. 3 we analyze the predictions of model Vj. In Sec. 4 we provide the new model Vn 
and analyze its predictions for the quark mixing angles and CP violation. Sec. 5 discusses 
a simple extension of model Vn to the neutrino and charged lepton sector and the resulting 
predictions. In Sec. 6 we have our concluding remarks. 

2 CP invariant SUSY Qq model 

2.1 Qq group theory and the Yukawa sector of the model 

We work within the context of supersymmetric standard model, with a non-Abelian fiavor 
symmetry acting on the three families of quarks, leptons and their superpartners. The 
group theory of is discussed in detail in Ref. [8]. We briefiy recall its salient features 
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relevant for model building. Qg is a binary dihedral group, a subgroup of SU{2), of order 12. 
It has the presentation 

{A, B; = E,B'^ = A\ B'^AB = A'^} . (2.1) 

The irreducible representations of Qq fall into 2, 2', 1, l', l", 1 ", where the 2 is complex-valued 
but pseudoreal, while the 2' is real valued. The {1, 1', 1", 1"'} singlets form a subgroup 
with the 1 and l' being real and the l" and 1 " being complex conjugates of each other. The 
group multiplication rules are given as 

1' X 1' = 1, 1" X 1" = 1', 1"' X 1"' = 1', 1" X 1"' = 1, I'x I'" = 1", I' X I" = I'" (2.2) 

2 X 1' = 2, 2 X 1" = 2', 2 X I'" = 2', 2' x 1' = 2', 2' x 1" = 2, 2' x 1"' = 2 (2.3) 

2 X 2 = 1 + 1' + 2', 2' X 2' = 1 + 1' + 2', 2 x 2' = 1" + 1"' + 2 (2.4) 

The Clebsch-Gordon coefficients for these multiplication can be found in Ref. [5]. 

In Table 1 we list the Qq assignment of the quark, lepton and Higgs chiral supermulti- 
plets in our modeljfl where Q,Q3, L, L3 stand for the SU{2)l quark and lepton fields, and 
if", H^, H'^, are the Higgs doublets. The SU{2)l singlet supermultiplets for quarks, 
charged leptons and neutrinos are denoted by m^, Wg, (i^, (ig, e'^, Cg and z^^, f|. Three pairs of 
Higgs doublets are introduced in order to generate fermion masses directly in the presence of 
Qq symmetry using renromalizable couplings. The singlet field T3 is needed to generate the 
Majorana mass for z/3. The other singlet scalar fields are needed to achieve spontaneous break- 
ing of Qe symmetry as well as CP symmetry without giving rise to pseudo-Nambu-Goldstone 
bosons. This point will be clarified in the next subsection. Table 1 also shows a fiavor uni- 
versal Z4 symmetry, the purpose of which is to realize an unbroken interchange symmetry in 
the scalar sector even after spontaneous symmetry breaking. Such an interchange symmetry, 
for which we have two solutions, enables us to predict one quark mixing parameter. 

Table 1: Particle content of the Qq model along with their transformation under QqX Z^^. 
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The most general Yukawa superpotential involving the quark and lepton fields invariant 
under the Qq x Z4 symmetry, assuming matter parity in the usual way, is: 

''Essentially the same model can be realized with any Q2N if N is odd and a multiple of 3. 
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+ {aeLselH^ + be{L * H'^)el + b'MH'^ * e") + Ce{L * e'')H^ + e ^ u} 
+ ^-^-^ + ^-3-3%, (2.5) 
where we have defined 

x-y = xiyi + X2y2 , X * y = xiy2 + X2yi , X -k y = xiy2 - X2yi . (2.6) 

We have used the exphcit basis for Qq given in Ref. and the notation u'^ = {—ul, U2) etc, 
for the right-handed Qq doublet fermion fields. Note that the Z4 symmetry plays no role in 
the construction of Eq. (12.51) . 

2.2 The Higgs sector 

In order to break the Qq symmetry spontaneously while avoiding pseudo-Nambu-Goldstone 
bosons one needs to introduce SM singlet Higgs fields. The minimal such set will involve a 2, 
2', 1' and two I's of Qq. These are listed in Table 1. The SM singlet S"s are needed to mix 
the Qq doublets iJ"''' with the singlets H^''^. Without the Qq doublet T there will be an 
accidental 0(2) symmetry in the Higgs potential. The 0(2) symmetry is violated by the cubic 
coupling of T. The field T3 is introduced for the Majorana mass for z/|, and the singlet 
U is introduced to generate a spontaneous CP violation and also to enable the spontaneous 
breaking of Qq x Z4 within the SM singlet sector. Thus the SM singlet Higgs sector employed 
appears to be the minimal set consistent with the demands we wish to meet. 

The most general Higgs superpotential involving the Higgs fields of Table 1 invariant under 
the Qq X Z4 symmetry along with the usual matter parity (with all the Higgs fields being 
even) has the form 

W^Higgs = Wu + Wst + Wh , (2.7) 

where 

Wu = fiuU^ + ^U^+ (Ai Si + A2 r| + AsT • T) U, (2.8) 
WsT = fiss Sl + fiTT-T + fiT, Tl + A'3 T ■ (T ® T) 

+A;[ - 2S2S{T^ + {Si - Sl)T2 ] + \'2S ■ ST; , (2.9) 
Wh = A'/ H'^{H^ * S) + A'2' (//" * S)H;^ + A'3' (i/" ■ H^)Ss (2.10) 

with the notation 

A -{3(^0) = A^{-B^Ci + B2C2) + A2iBiC2 + B2C1) . (2.11) 

Thus T-{T®T) = STiTl -Tf. The Z4 symmetry has restricted the form of Eqs. fl^ - fl^TTU]) : 
without the Z4, the following couplings would be allowed: 

<iggs = (-i^r^i + H:^H^)T, + {H^Ht + H:^Ht)T2 . (2. 12) 
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We wish to avoid these terms, since in their absence we can define an unbroken discrete 
symmetry, as discussed below. 

The Higgs potential contains F terms derived from Eqs. 02.81) - f|2.10p . D terms associated 
with SU(2)l X U{1)y breaking, and the following soft SUSY breaking Lagrangiaij^l 

Csoit = m^|f/p + m|(|5i|2 + |52|')+m|3|^3|' + m|(|Ti|2 + |T2p) + m|3|T3|2 

+ {BuU^ + Bs, SI + BtT-T + Bt, 

+ [AU^ + A, Si + A2 T| + A3 (T ■ T)] U + A'^T-{T^T) 

+ A[[ - 2S2SiTi + {Si - Sl)T2 ] + A^S ■ ST^ 

+ A[ H'^{H^ * S) + A'^ (if" * S)H^ + A'^ {H^ ■ H^)S-i + h.c. ] , (2.13) 

where the ■ and * products are defined in (12.61) . We assume CP invariance, which implies 
that all the Yukawa couplings and the parameters in the Higgs potential are real. The 
Higgs potential would then admit two interesting minima which leave two separate discrete 
symmetries Vi or Vu unbroken. We analyze these two ground states in the next two sections. 



3 Ground state with unbroken interchange symmetry 

The following symmetry Vi is respected by the Qq x Z4 invariant Higgs superpotentials Eqs. 
dZSD - fl210l) . and fl27[3|l and the D terms: 

HI ^ HI H^ ^ Hi Si ^ S2, T2 ^ -T2, 
H^ ^ HI Hi -> Hi S3^Ss,Ti^ Ti, T3 -> T3, U. (3.14) 

The VEVs of the various Higgs fields can be consistently chosen such that this symmetry 
remains unbroken: 

^ii^'^) = (^Hl""^ = vl'''e"i'+\ (hI'^'j = vl''e'^^^'\ {Si) = {S2) = vse"^^ , 
(Ti) = vre^'^^ (T2) = 0, (^3) = vs.e^'^'^, (T3) = vr.e'^^^, {U) = we'^^ . (3.15) 

In Eq. (13.151) . we have explicitly displayed the complex phases. It should be noted that this 
symmetry Vi is an accidental symmetry of the Higgs potential, and is not respected by the 
full theory. For example, the Yukawa sector explicitly breaks this symmetry. Nevertheless, 
the existence of Vj enables us to choose a ground state given as in Eq. (13.151) consistently. 

^We have used the same symbol for the scalar components as the superfields. 
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We have explicitly verified that the minimum of Eq. f l3.15p is indeed a local minimum, 
and that spontaneous breaking of Qq x Z4 and CP symmetries occurs without generating 
pseudo-Nambu-Goldstone bosons. The scalar spectrum of our model is in fact arrived at by 
meeting these requirements. 

In the ground state Vi, the mass matrices for the up and down quarks take the form: 



/ ^e^A0„,,x 



u,d 



V2 



V2 ^ V2 

Here we have defined the following parameters: 



0/ rrf 



(3.16) 



Au,d = au,d v^'"^, Bu4 = V2 bu,d f r'^, B'^,d = K,d ^i''^; Cu,d = Cu,d v^''^, 

A0.,d = 03'' - 01-' • (3.17) 

We have ignored irrelevant overall phases of the two mass matrices. CP invariance of the 
Lagrangian implies that the parameters (A„ rf, B^^, B'^ ^, Cu,d) are all real. In this case, 
after a common 45 degree rotation in the (1-2) sector that would set the (1,3) and (3,1) 
entries of Mu,d of Eq. f l3.16p to zero without inducing CKM mixing, we can write 

= Pu,dMu,dPu4 , (3.18) 

where Mu^d are real matrices given as 




Mu,d= I -Cu,d Bu,d I , (3.19) 

and Pu,d are diagonal phase matrices given as 

Pu,d = diag.{e-*^'^-■^ e*^'^"■^ 1} . (3.20) 
The CKM matrix is then given by 

VcKM = OlPOd , (3.21) 
where Ou,d are the orthogonal matrices that diagonalize M^^d via 

OldMu,dMlaOu,d = diag.{m2_^, m^^, mjj , (3.22) 
and P is a diagonal phase matrix 

P = diag.{e^'^, e-'^, 1} (3.23) 
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with (f) = A(f)d — A(j)u. Since Mu and each has four real parameters, once the six quark 
masses are fixed, 0„ and Od will have one undetermined parameter each. These two param- 
eters and the phase appearing in the matrix P of Eq. f l3.23p will completely fix the three 
CKM mixing angles and the one CP violating phase. That will lead to one sum rule involving 
the CKM mixing angles, the CP violating phase, and the quark mass ratios. This prediction 
was analyzed in Refs. [8l[T2] and shown to be fully consistent with data. 

Here we update the results of Ref. |S] for the quark mixing parameter prediction. We use 
the most recent values of the quark masses. Lattice calculations have reduced the errors in 
the light quark masses, which we adopt for our fits. Furthermore, we compare the model 
prediction with the global fits provided recently in Refs. [9],[T3] assuming specific new physics 
contributions. The new physics contributions are motivated by certain discrepancies that 
have been observed in the CKM fits. For example, the CP violation parameter r/ determined 
from sk differs from that obtained from the decay B — > J/ipKs by more than 2 standard 
deviations. We compare our model fits with the best fit of the standard model, as well as 
with the best fit for Scenario 1 of Refs. [HllIS]- This scenario is characterized by independent 
new contributions to B^^s — B^.s mixing amplitude. It turns out that there is room for 
small (~ 25%) new contributions to these mixings in our model, arising from gluino-squark 
box diagrams. The Qq assignment of quarks implies that the third family squark is not 
degenerate with the first two family squarks (which are nearly degenerate). Once the quark 
mass matrices are diagonalized, there will be small off-diagonal entries in the squark mass 
matrix, which leads to Bd,s — Bd,s mixings. These diagrams have been evaluated in Ref. [TT] . 
While real, these amplitudes are still in the interesting range for new physics to influence the 
CKM parameter flts. In Ref. [H] the radiative corrections to these mixing parameters, arising 
through Higgs boson exchange, have been computed, and have been shown to be complex. 
Thus, it appears that the Qq model admits small deviations in the CKM flts to Bd^s ~ Bd^s 
mixings. It should be noted, however, that the prediction of the present model agrees well 
with the best flt values of the CKM flts, with or without new physics assumed. 

Guided by the analytic expressions for the CKM mixing parameters from (13.191) - (13. 23 p . 
we have done a numerical flt to all quark masses and mixings. An excellent flt is obtained 
with the following choice of parameters at = 1 TeV: 

Au/mt = 0.9963, Bjmt = 0.06051, Bjmt = 0.06051, Cu/mt = 1.748 x 10-^ 
Ad/rub = 0.8895, Bd/rub = 0.04214, B'Jnib = 0.4554, Cd/rub = -5.043 x 10'^, 

(p = 0.71875. (3.24) 

The resulting mass eigenvalues at /i = 1 TeV are: 

m„ = 1.25 MeV, = 552 MeV, 

rud = 2.74 MeV, = 50.0 MeV, (3.25) 
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where we have used rrit = 150.3 GeV and mb = 2.46 GeV. These values are to be compared 
with quark masses extrapolated from low energy scale to /i = 1 TeV 



= 0.85 ~ 1.65 MeV , = 2.05 ~ 2.90 MeV , 
ms = 39.6 ~ 64.4 MeV , = 502 ~ 570 MeV , 
mb = 2.39 ~ 2.53 GeV , rut = 148.9 ~ 151.6 GeV , (3.26) 

where we have updated the result of [15j by using the updated quark masses given in PDG 
2011 |il6j, while neglecting the uncertainties due to the RG running. The input values of Eq. 
fl3.24p give also the following output for the CKM parameters: 

A = 0.2252, A = 0.7962 , p = 0.1613 , 77 = 0.4230 , 
sin2/3 = 0.8042 , a = 84.1 [deg] , /3 = 26.8 [deg] , 7 = 69.1[deg] , (3.27) 

which should be compared with the fit result of the CKMfitter group (scenario I) [9] 

A = 0.22542 ± 0.00077 , A = 0.801^[J:[J?^ , (3.28) 

-p = 0.159^^1 , V = 0.438in^^ , sin2/3 = 0.813^^^^ , 

a = 79tll [deg] , P = 27.2+^:^ [deg] , 7 = 70.0+^5 M • (3.29) 

Since there are nine model parameters for six quark masses and four CKM mixing param- 
eters, we can make one prediction in a two dimensional plane if we fix eight of the nine model 
parameters. To fix these eight parameters we use the quark masses, A and A given in f l3.26p 
and f l3.28p . respectively. Fig. [1] shows the prediction in the p — f] plane, and Fig. |2] shows the 
prediction in the /3 — 7 plane. The CKMfitter group best fit values fl3.29p are also indicated 
in these figures. We see from Eqs. f l3.25p . Fig. [T]and Fig. [2] that the model Vj reproduces the 
quark masses, CKM mixings and the CP violating phase in an excellent way. 

4 New ground state with unbroken interchange sym- 
metry Vii 

The same Higgs potential as derived from Eqs. (12. 7p . and the soft SUSY breaking Lagrangian 
fl2.13p including the D terms, with all parameters taken to be real so that CP is an exact 
symmetry, admits a new unbroken interchange symmetry as given below: 

H^*, S3 -> S;, Ti ^ T*, T^^T;, U ^U*. (4.30) 
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Figure 1: The prediction in the p — fj plane for the model Vj, where we have used as the input parameters; 
the quark masses, A and A given in Eqs. (j3.26p and p.28p . respectively. We also have imposed the constraints 
on the quark mass ratios |16) : 2171^/(171^ + md) = 22 ~ 30, nis/md = 17 ^ 22, mu/md = 0.35 ~ 0.60. The 
crosses are the CKMfitter group values [S]; blue (scenario I) and red (SM). 



This symmetry Vu enables us to choose a ground state given by 
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(4.31) 



where the complex phases are all explicitly displayed. Note that there are only three phases, 
(ps, (pu and (pd in the VEVs, along with a purely imaginary VEV of T2. 

In the background Vu, the fermion mass matrices M^^d following from Eq. (12.51) take the 
form 
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~Cu.d 



'^u,d .Pi C \ 



B' 

V2 ^ 





B' 

V2 ^ 



V2 
^/2 ^ 



(4.32) 



with the parameters as defined in Eq. f l3.17p . CP invariance of the Lagrangian implies that 
the parameters {A^^, Bu,d, B'^^^, Cu,d} are all real. 
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Figure 2: The prediction in the /? — 7 plane for the model Vi- The input parameters and the constraints 
are the same as for Fig. [1] 

Model P//, while different from model "P/, is just as predictive in the quark sector as Vi. 
It is then interesting to see if the quark mixing sum rule of Vn is consistent with data. To 
address this question we proceed to diagonalize M^^ of Eq. f l4.32p . The phases in the matrices 
of Eq. fl4.32p can be factorized: 

M,,rf = P„,rfM;,P,,d (4.33) 

where 

Pu4 = diag.{e*<^-^ e-'^"•^ 1} (4.34) 

with M^' ^ given as in Eq. (14. 32 p . but with (j)u,d set to zero. Quark field redefinitions can absorb 
the phases in Pu^d, however a phase matrix will then appear in the quark mixing matrix: 

P = diag.{e'^, e"^^, 1} , (4.35) 

where 

<P = <Pd-<Pu. (4.36) 

Now we do a 45 degrees rotation in the (1-2) plane to bring M^^ into M„,d as given in Eq. 
fl3.19p . but this will generate a non-trivial quark mixing matrix given by 

(cos (j) i sin (f) \ 
isin0 COS0 . (4.37) 
1/ 
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The CKM mixing matrix is then obtained as 

^CKM = O^KOd , (4.38) 

where Ou,d diagonahze the matrices of Eq. fl3.19p as specified in Eq. f l3.22p . 

Using the approximate analytic expressions for the CKM mixing parameters, we have done 
a numerical fit to all quark masses and mixings within this model. An excellent fit is obtained 
with the following choice of parameters at /x = 1 TeV: 

Au/mt = 0.01389, B^/mt = -0.003282, B'Jrrit = 0.9999, Cu/mt = 1.381 x 10'^, 
Ad/nib = 0.9020, Bd/nib = 0.04512, B'Jrrib = 0.4297, Cd/nib = 4.554 x 10~^, 

= 0.1038. (4.39) 

The resulting mass eigenvalues at = 1 TeV are: 

m„ = 1.12 MeV, = 535 MeV, 

ma = 2.27 MeV, m, = 50.0 MeV, (4.40) 

where we have used rrit = 150.3 GeV and mb = 2.46 GeV as in the case of Vi. These values 
are to be compared with quark masses given in Eq. fl3.26p . The input values of Eq. f l4.39p 
give the output for the CKM parameters: 

A = 0.2254, , A = 0.7987 , p = 0.1575 , r/ = 0.4231 , 
sin 2/3 = 0.8021 , a = 83.7 [deg] , (3 = 26.7 [deg] , 7 = 69.9[deg] , (4.41) 

which should be compared with the fit result of the CKMfitter group (I3.29p . 

Fig. 13] shows the prediction in the p — f] plane, and Fig. S] shows the prediction in the 
/3 — 7 plane for model Vn. The CKMfitter group best fit values fl3.29p as well as the SM best 
fit values are indicated in these plots. As in the case of Vi, we see from Eqs. (14. 40 p . Fig. [3] 
and Fig. H] that the model Vn also reproduces the quark masses, CKM mixings and the CP 
violating phase in an excellent way. 

5 Predictive scenario for neutrino mixing 

The lepton sector of model Vi with the Qq assignment given in Table 1 has been studied 
in Ref. [8], and therefore we will not discuss it further here. It is interesting to see if there 
are any constraints on neutrino oscillation parameters for model Vn. Here we explore an 
alternative possibility of the Qq assignment for the leptons, which is given in Table 2. 
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Figure 3: The prediction in the p — f] plane for the model Vu, where we have used as the input parameters; 
the quark masses, A and A given in Eqs. (j3.26p and p.28p . respectively. We also have imposed the constraints 
on the quark mass ratios |16) : 2ms/{mu + md) = 22 ~ 30, ms/md = 17 ~ 22, m„/md = 0.35 ~ 0.60. The 
crosses are the CKMfitter group values [3; blue (scenario I) and red (SM). 

Table 2: An alternative x assignment for the leptons. 
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In this new scenario, the leptonic part of the superpotential fl2.5p becomes 

+ ^v^-v^^'^-^vlvlT,, (5.42) 

where the ■ and ® products are defined in (12.61) and (12.111) . respectively. 
The Majorana mass matrix for the right-handed neutrinos is given by 



Ml 











Ml 











M3 



M^c = Ml , (5.43) 
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Figure 4: The prediction in the /? — 7 plane for the model Vu. The input parameters and the constraints 
are the same as for Fig. [3] 

where M3 = a^c vt^. Note that Mi and M3 are both reaL The Dirac neutrino and charged 
lepton mass matrices are: 

M^D = I ae-*"^" C^e^'f'- j , = [ Cie'''^'' Cte''^^ Bte'"'^^ j .(5.44) 
\ B^'^-- 5^e-^*" aJ \ B'^e'"^" B'^e-'"^" / 

The hght neutrino Majorana mass matrix is found (by the seesaw formula) to be 

(2p|cos(20„) -2ip2P4 sin(20„) \ 

2picos(2(/.J 2p2P4 , (5.45) 
-2ip2P4: sin(20„) 2p2P4 -pj + "^pI cos(20„) / 

where 

pI = {C^f/Mi , pI = -{A^f/M, , pI = {Blf/M, . (5.46) 

We have assumed that Mi is positive, while M3 is negative. When 0„ = 0, the neutrino 
mass matrix is exactly the same as the matrix discussed in [T7], and yields only a tiny 
Ue3 ~ nie/m^ ~ 10~^, where f/gs is the (e, 3) element of the Maki-Nakagawa-Sakata (MNS) 
mixing matrix [/mns- It was also shown there that the mass matrix (15.451) with (pu = ^ can 
yield consistent neutrino masses and mixing only if M3 is negative, and the mass spectrum is 
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inverted. This conclusion also applies to the present case with non- vanishing For non-zero 
(pui we obtain U^z ~ sin2(/)„, which can be small or large. We vary \U(.z\ in its entire range 
allowed by experiments and correlate its value with other observables. 
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Figure 5: The prediction in the |C/e3p — i^mns plane for the model Vn with the Qq assignment of the leptons 
given in Table 2, where we have used the parameters given in (|5.52p . and (j>d — (f'u + 0.1038. The dashed 
vertical line is the maximal CP violation. 

We make the matrix fl5.45p real by redefining ui = w^. The resulting mass matrix M]^^"^^ 
can be diagonalized by an orthogonal matrix Oy as O'^M^^^^^O^. As for the charged lepton 



mass matrix Mi, we can obtain hierarchical masses, e.g., rrie ~ 5^ , 
Keeping this in mind we rotate according to 



rrir 



M, 



PlM^Pr 



where 








R 




~«<Pd 



(5.47) 



(5.48) 











t . 



Then we consider M^M^ in the limit B[ — t- 0, (i.e. nie — )■ 0), and find 

C|(3 - cos(40,)) C|sin(40rf) 

C|sin(40rf) 25| + C|(l + cos(40rf) 

0, 



MpM, 



t 



(5.49) 
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Figure 6: The prediction in the |C/e3p— < ^ee > plane for the same input parameters as Fig. [5] 



The eigenvalues in the hmit are 



, ~ C|(3 - cos(40rf)) , ml - IB] + + cos(40rf) 



(5.50) 



and the (inverse) diagonahzing orthogonal matrix (OjMiMi Oi) is found to be 





1 



\ ^ sin(40rf) 



-^sin(40d) 



(5.51) 



Since the relative phase = 0^ — 0u is fixed in the quark sector, there are seven independent 
parameters in the lepton sector. We use [18] : 



rrie 



0.511 MeV , rrif, = 105.7 MeV , = 1.777 GeV , \U, 



e2 



0.318 



+0.019 
-0.016 



Am?3 



(2.40 t°:}?) X 10-=^ eV^ , Ami, = (7-59 t^ofs) x 10"^ eV^ 
as input parameters. The MNS neutrino mixing matrix is then given by 

f/MNS = OJPlPuO, X diag.{ 1,^,1}, 



(5.52) 



(5.53) 



where the last phase factor multiplied with Uum is the Majorana phase, and Pi, = diag.{ i , 1 , 1 
which was introduced to make the matrix fl5.45p real. In the lepton sector we have only one 
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^— 0.54 
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0.46 
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e3 

Figure 7: |C/^3p against |C/e3p for the same input parameters as Fig. [5] 
free phase (pu, which controls f/es- In the following calculations we use (pd = 4'u + 0.1038 (see 

Fig. [5] shows the Dirac phase (5mns (in the convention of Ref. [H]) against |t/e3p. We see 
that the model predicts nearly maximal CP violation. This can be understood as follows. 
Consider the limit me, 0^1, 0^ — )■ 0. In this limit, only P^, contributes to Smns, and the first 
element of P^, e*'"/^, appears as the Dirac phase. 

It is possible to predict the effective neutrino mass < niee >= l^u^U^i + "^1/2^^62 + "^!^3f^e3l 
for neutrinoless double beta decay as a function of \Ue3\. Note that the first row of OJPlPu 
is diag.{ , , 1 } in the mg — )■ limit. Since Ou is real, the first and third elements of the 
first row of Uuns are real, while the second element is purely imaginary. Therefore, 

< > — Irriiy-^ cos^ ^soi — nT'U2 sin^ ^^soi| — ^^1^2 cos20soi — 0.4 m^^. (5.54) 

In Fig. O we plot the prediction in the |f/e3p— < nT-ee > plane, which verifies the rough 
estimate above. The main contribution to lU^^l comes from Oe. In the limit me,(f)u 
0, it is exactly 1/a/2, so the maximal mixing. The deviation from the maximal mixing 
has terms proportional to rrie/m^ and to sin20„. In Fig. [7] we plot |t/^3p against |t/e3p, 
verifying our expectation. Note that the entire range of {Uesl allowed by experiments currently 
is also allowed by atmospheric neutrino oscillations. But once the {Uesl is measured, the 
model will make precise prediction for |f/^3| which can be scrutinized with improved precision 
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experiments. 

6 Conclusions 

The Qq model of flavor is constructed to solve the SUSY flavor problem of the supersymmetric 
standard model. It also yields an interesting prediction for the quark mixing parameters, 
which compares very well with experimental data. An unbroken interchange symmetry plays 
an important role in obtaining the quark mixing parameter prediction. In this paper we have 
updated this prediction, and compared it with the best flt values within the standard model 
as well as with new physics contributions assumed in B^ ., — B^^s niixing amplitudes. The 
model prediction is in very good agreement with the data. 

A major observation of the present paper is the existence of a new minimum that violates 
CP symmetry spontaneously, but leaves a new interchange symmetry unbroken. In this 
minimum, there is again a prediction for quark mixing parameters. We have analyzed this 
prediction and found that it flts data (within the CKM model and with new physics included) 
rather well. We have extended this symmetry to the leptonic sector, and have found various 
correlations between neutrino oscillation parameters. 

We conclude with several comments on the new solution found. 

(1) The SUSY flavor problem is solved in the new ground state Vu in the same way it is 
solved in Vj. Qq invariance requires the flrst two family squarks and sleptons to be degenerate 
in mass, which provides the needed SUSY GIM mechanism. Since after Qq breaking the Qq 
doublet and singlet quark states mix, there is residual flavor violation mediated by the SUSY 
particles, but such FCNC processes are well within experimental limits. 

(2) The SUSY CP problem is solved in the model by virtue of spontaneous CP viola- 
tion. The fundamental parameters in the Lagrangian are all real, complex phases develop 
only spontaneously via the VEVs of if"'*^ and S, T, U flelds. This implies that the soft SUSY 
breaking parameters such as the gluino mass are all real, which alleviates bulk of the SUSY 
phase problem. The trilinear SUSY breaking A-terms are not proportional to the correspond- 
ing Yukawa couplings, however the phases in these A-terms, since they arise spontaneously, 
will align with the phases in the fermion mass matrices. Thus the A-terms do not generate 
CP violation. There is CP violation arising from the /i-terms, but as suggested in Ref. [TT], 
if the Higgsino masses are parametrically smaller than the squark and slepton masses, this 
CP violation is not excessive. We also note that in the new minimum Vn-, the spontaneously 
induced phase that is necessary for KM CP violation is rather small, ~ 0.1038. One can then 
assume an approximate CP symmetry for the entire Lagrangian, where all the phases remain 
small, of this order. This will further suppress the SUSY phase effects. 

(3) The new interchange symmetry Vu might appear to be CP transformation, but it 
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actually is not. If it were CP transformation, when extended to the fermion Yukawa sector, 
that would make the parameters Cu^d/,u in Eq- (12. 5 p purely imaginary. CP violation will 
then disappear from the CKM matrix, as it should, since this symmetry remains unbroken. 
The symmetry Vn is an accidental symmetry of the Higgs potential, and is not respected 
by the Yukawa couplings, just as it was for the interchange symmetry Vi. This state leads 
to a new sum rule involving the quark masses and CKM mixing parameters, which is found 
to be in good agreement with data. Extension of the model to the neutrino sector, by 
changing the assignment of the leptons, can lead to a predictive scenario. In this version 
we find that neutrino mass hierarchy is inverted with nearly maximal CP violation along with 
nearly maximal mixing of atmospheric neutrinos. Thus the model lends itself to experimental 
scrutiny in the near future. 

(4) The question of whether it is possible to obtain a large CP violation in the — 
mixing for the case of Vn, as in the case of Vi [H], remains to be studied. To distinguish 
two ground states of the same model, precise measurements of the CKM parameters [19] and 
neutrino oscillation parameters |2U] as well as precise determination of the quark masses [2T] 
are indispensable. 
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